home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / TTEST.DEM < prev    next >
Text File  |  1991-04-29  |  1KB  |  50 lines
  1. PROGRAM d13r3(input,output);
  2. (* driver for routine TTEST *)
  3. (* generate gaussian distributed data *)
  4. CONST
  5.    npts=1024;   (* npts must be > mpts for TTEST.PAS *)
  6.    mpts=512;
  7.    eps=0.03;
  8.    nshft=10;
  9. TYPE
  10.    glnarray = ARRAY [1..npts] OF real;
  11.    glmarray = ARRAY [1..mpts] OF real;
  12.    glnparray = glnarray;
  13. VAR
  14.    glinext,glinextp : integer;
  15.    glma : ARRAY [1..55] OF real;
  16.    gliset : integer;
  17.    glgset : real;
  18.    data1 : glnarray;
  19.    data2 : glmarray;
  20.    i,idum,j : integer;
  21.    prob,t : real;
  22.  
  23. (*$I MODFILE.PAS *)
  24. (*$I AVEVAR.PAS *)
  25.  
  26. (*$I GAMMLN.PAS *)
  27.  
  28. (*$I BETACF.PAS *)
  29.  
  30. (*$I BETAI.PAS *)
  31.  
  32. (*$I RAN3.PAS *)
  33.  
  34. (*$I GASDEV.PAS *)
  35.  
  36. (*$I TTEST.PAS *)
  37.  
  38. BEGIN
  39.    gliset := 0;
  40.    idum := -11;
  41.    FOR i := 1 to npts DO data1[i] := gasdev(idum);
  42.    FOR i := 1 to mpts DO data2[i] := (nshft DIV 2)*eps+gasdev(idum);
  43.    writeln('shift':6,'t':8,'probability':16);
  44.    FOR i := 1 to (nshft+1) DO BEGIN
  45.       ttest(data1,npts,data2,mpts,t,prob);
  46.       writeln((i-1)*eps:6:2,t:10:2,prob:10:2);
  47.       FOR j := 1 to npts DO data1[j] := data1[j]+eps
  48.    END
  49. END.
  50.